/ y 

V ' - 

Final Technical Report / y ~ 7 

Aeroelastie Stability and Response of 

Rotating Structures 

NASA Grant NSG-3139 


Theo G. Keith, Jr. 

epartment of Mechanical Engineering 
University of Toledo 

N93-16560 
Unc 1 as 
G3/08 0140846 

January 1993 


CNASA-CR-191803) AEROFIacttc 

^ R0TAT ^G 

(Toledo Univ.y 7 p C Rep ° rt 



Aeroelastic Stability and Response of 
Rotating Structures 


Summary of Grant NSG-3139 

A summary of the work performed under NASA grant number NSG-3139 is presented 

S/r d f tailS °‘ th f 5 efforts can in the cited references. Thfs grant has 

1 d i J h Z' develo P ment of analysis methods for predicting loads and instabilities of 
wmd turbmes and the initiation and development of three new areas of “h o a i d 
the Advanced Turboprop Project (ATP). u 


(A) Analysis of Wind Turbines 

This work was started to assist the Department of Energy (DOE) to analyze and desien 
efficient wind turbines for extracting wind energy. Initially, two versions of § a 
commeraally developed computer program (MOSTAS) for predicting loads and 

bladpH ^ ° wmd turbmes, were evaluated by comparing with measured data for two 
bladed machines. It was concluded that the version having the periodic coefficients 
solution performed better than the version that used time averaging. Refs. A.l and A 2 
At the same time, a simple model with five degrees of freedom, was developed for 
hl^H 8 ; T 6 Whir ‘ ! Utter stabili ‘y of horizonfal-axis wind turbines wiZa two- 

b d d ! et T ng r ° t0r ' ThlS model accoun ted for the out of plane bending motion of 
each blade, the teetering motion of the rotor, and both the pitching and yawing motions 
of the rotor support. The results showed that the DOE/NASA Mod-2 design^vas free 
from whirl flutter. Ref. A.3. A parametric study was performed showing the Effect of 
variation 8 in rotor support damping, rotor support stiffness, and the effecf of pitch-flap 
coupling on pitching, yawing, teetering, and blade bending motions. P P 

SeC f °? d degree non-linear governing equations for the analysis 
vertical axis and horizontal axis wind turbines were developed. Refs A 4 and A 5 

^n? US6d t0 inV6Stigate ^ **** of many practical 

flexible blades beCaUSe ° f lar S e deformations of the 

assumDtinn nf ^ ordering scheme, which is consistent with the 

assumption of a slender beam, was used to discard some higher order elastic and 

accountedTo/bo^h Seco ^ d : degree nonlinear equations. The blade aerodynamic loading 

based nn * n b tb ! nd Shear and tower shadow ; it was obtained from strip theory 

afffofl Teor^ THp 6 ir PPr ° Ximati0n ° f tw °- dime ™°" a l incompressible, unsteady 
airtoil theory. The resulting equations had periodic coefficients. y 

vibration and th j rd degree geometric nonlinear terms on the 

W a S studied in l bl dy oharactenstics of rotating, pre-twisted, and pre-coned blades, 

mPthnd Tb ? efS ‘ A-6 and A ' 7 ' The et l uatlons were solved using the Galerkin 
and aim 1 Ctl 0 j S were compared with those obtained from MSC / NASTRAN 

nd also with measured data. It was shown that the spurious instabilities, observed for 



thin, rotating blades when second degree geometric nonlinearities were used, could be 
eliminated by including the third degree elastic nonlinear terms. The inclusion of third 
degree terms improved the correlation between theory and experiment. A study to 
investigate the vibration and buckling of rotating, pre-twisted, pre-coned beams 
including Coriolis effects was performed in Refs. A.8-A.10. It was shown that the 
oriolis effects were necessary for blades of moderate-to-large thickness ratios. Thus, 
the linear Coriolis terms associated with pre-cone could be neglected in the dynamic 
analysis of advanced turboprop (propfan) blades. The results also showed the 
possibility of buckling due to centrifugal softening terms for large values of pre-cone 
and rotation. During this study, an improved finite difference method was also 
developed; this improvement eliminated most of the shortcomings associated with 
conventional first order methods and also provided faster convergence, Refs. A.11-A.13. 

(B) Aeroelastic Analysis methods for cascades including blade and disk flexibility 

The first new area of research aimed at helping the Advanced Turboprop Project is 
cascade aeroelasticity i.e., an investigation of the response and stability of multibladed 
structures, such as propfans, turbines, and compressors. This is an important topic in 
the design of propulsive elements of any engine. The work started with the pioneering 
paper on the development of methods for the stability analysis of mistuned rotors. 
Ref. B.l. This paper used only a typical section structural model, with bending and' 
torsion motion, in incompressible flow. However, this formed the basis for all the 
methods developed to date. The aerodynamic models were later improved to include 
the effects of compressible flow in general and supersonic flow, in particular. Refs. B.2 
and B.3. Even this simple model produced considerable insight into stability and 
response aspects of cascades. The study showed that mistuning is beneficial for the 
stability of the cascade, but not for its response. Since mistuning exists naturally, this 
has to be studied very carefully. The study resulted in the development of computer 
programs TUNEBT and MISER, which are widely used even today in the design and 
testing of the stability and forced response of cascades. 

The later studies concentrated mainly on the development of the structural model. The 
typical section structural model was replaced by a beam model. Refs. B.4-B.8, and a 
finite element model. Ref. B.9. The available two-dimensional theories were used in a 
stripwise manner. The structural model was extended to include disk flexibility, which 
required a special formulation to include the effect of both rotating and non-rotating 
parts. Ref. B.10. The corresponding computer programs, BEAM, ASTROP3, and BEDE 
are now routinely used for further research. During this period, methods to automate 
flutter calculations were also developed. These methods drastically reduced the time 
required to calculate the flutter conditions of a cascade, propfan or of a turbomachine. 
Ref. B.ll. All the developments until 1988 were included in the survey paper presented 
at Lewis Structures Technology, Ref. B. 12 . 



(C) Stall flutter analysis 

»™l a bn it y In repitedtotleTonf tZS'oI pr ° pfanS ' 

SES^^ 

ONESwSrr Were reviewed ' Ref - CI 11 was concluded that the 

of Qfaii ac £ aI y ses ' Refs. C.2 and C.3. It was found that the models predicted the onset 
a,rfnnl' ° bserved “ experiments. However, static stall data for NACA 16 senes 

ACTROP2 s T II q T U,r h T g0 °5* correlation - The associated computer proglam 
m QP2 ‘ S I ALL was dellvered to Hamilton Standard for evaluation and use* The 

03 WaS awarded the S “* Branch b 0 est a p n ape U r Se of*: 


(D) Computational Aeroelasticity 

The third area of research, which was started under this grant and led to leading pHctp 
technology, is computational aeroelasticity. It was known that linp^r a & ^ 

Rories used in the cascade aeroelastic studies ^LTa^nZ 

effects of airfoil shape angle of attack, and transonic flow conditions Only through 
computational methods, could these effects be accounted for properly In^tiallv m 
become familiar with the broad area of computational aeroeLtirity a computer 
program was developed to study isolated airfoils. The objective of the^studv w^s to 

vfscositv thi S cknes/a W nd h '° inves,i 8 ate the <**•» °< roladonai flow,' 

metho1iTLl h c« CO H PUta I‘ 0naI * imeS and f ° ° btain a basic understanding of analysis 
methods for cascades, a time accurate full potential solver was develoned and „J ^ 

L7e 8 w a as U stf ^ ^ a tyP^Sn^rS 
data, a"ol«o^^^ 
methods are time consuming, e?ficL, 

TheIf U ‘ Tu tlme 1x1 r calculating the unsteady aerodynamic coefficients Ref D4 

d„ P m ametri H f d 7 ndUdin S comparison 'ones^Us 

thesis. Ref. d 2 ^ d ' d me domain method s, were presented in a Ph.D. 
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